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H , A new ansatz is presented for a Lax pair describing systems of parti- 

cles on the line interacting via (possibly nonsymmetric) pairwise forces. 
Particular cases of this yield the known Lax pairs for the Calogero-Moser 
and Toda systems, as well as their relativistic generalisations. The ansatz 
leads to a system of functional equations. Several new functional equa- 
tions are described and the general analytic solution to some of these is 
given. New integrable systems are described. 



Abstract 



1 Introduction 

Completely integrable systems arise in various diverse settings in both mathe- 
matics and physics and accordingly have been studied from many different points 
of view, a fact which underlies their importance and interest (see for example 



||23|| .) Within this area the study of Lax pairs, a zero curvature condition, plays 
an important role. The construction of such Lax pairs has followed many routes 
and this paper will further investigate the connection between functional equa- 
tions and such zero curvature conditions. The essential idea in this approach is 
to reduce the constraints of the Lax pair L, M implicit in L = [L, M] to that 
of a functional equation. Our study will broaden the ansatz for the Lax pair 
and correspondingly lead to a more general functional equation than hitherto 
studied. This enables us to understand many of the known integrable systems 
(and their corresponding functional equations) from a unified perspective. The 
symmetries of the functional equation we obtain are very large and this group 
relates distinct functional equations and their corresponding physical models. 
We feel this connection between functional equations and completely integrable 
systems is part of a broader and less well understood aspect of the subject. 
Functional equations have of course a long and interesting history in connection 
with mathematical physics and touch upon many branches of mathematics R |3| . 
Novikov's school for example considered the Hirzebruch genera associated with 
the index theorems of known elliptic operators and showed that these arose as 
solutions of functional equations. More recently Ochanine showed the string in- 
spired Witten index could be described by Hirzebruch's construction where the 
the functional equation was that appropriate to an elliptic function. These same 
functional equations arise (as we shall later see in more detail) in the context 
of completely integrable systems. The latter also appear in the study of con- 
formal and string theories and this connection between string theory and finite 
dimensional completely integrable systems needs to be better understood. The 
functional equations and integrable systems we shall discuss arise naturally in 
investigations of the KP and KdV equations [I], |2T[ . 

To make matters concrete let us consider how such functional equations arise 
in the context of integrable systems of particles on the line. Here one starts with 
an ansatz for the matrices L and M of the Lax pair and seeks the restrictions 
necessary to obtain equations of motion of some desired form. These restrictions 
typically involve the study of functional equations. The Calogero-Moser ]T5| 
system provides the paradigm for this approach. Beginning with the ansatz (for 
n x n matrices) 

L jk = Pj5jk+ 9 (l-5 jk )A(q j -q k ), M jk = g [8 jk y^ B(q j -qi)-(l-S jk )C(q j -q k )] 

1+3 



one finds L = [L, M] yields the equations of motion for the Hamiltonian system 
(n > 3) 

H = — y^Pj + g 2 2_. U(qj — Qk) U(x) = A(x)A(— x) + constant 

3 3<k 

provided C(x) = —A'(x), and that A(x) and B(x) satisfy the functional equation 

A(x)A'(y) - A{y)A\x) = A(x + y)[B(x) - B(y)]. (1) 

The solutions to this functional equation may be expressed in terms of elliptic 
functions and their degenerations.^ Krichever used this functional equation in 



his proof of the 'rigidity' property of elliptic genera [22] and it appears when 



discussing rational and pole solutions of the KP and KdV equations [21, M. Dif- 



ferent starting ansatz lead to the relativistic [28, H, 29] Calogero-Moser systems, 



the Toda J53, |23| and relativistic Toda S, [| equations. Underlying the corre- 



sponding functional equations of each of these models lies the addition formula 
for elliptic functions. Further, in studying the quantum mechanics of these sys- 
tems similar functional equations arise when factorising the ground state wave 



function 13 



In this paper we will introduce a new ansatz that includes the previous ex- 
amples (together with their functional equations) as special cases. We shall be 
seeking Lax pairs that lead to equations of the form 

q) = y^(q + bqj)(a + bq k )V jk (qj - q k ). (2) 

k+3 

Here we are allowing the interaction^ Vj k to in principle depend on the choice of 
pair j, k: when Vj k is the same for all pairs we have a system of Calogero-Moser 
type while if Vj k is the same for pairs j, j ± 1 and zero otherwise we have a Toda 
system (of A n type). When the constant b = we obtain the nonrelativistic 
systems and when a = we obtain the relativistic systems previously exam- 
ined. Indeed, because the interactions of the system only depend on coordinate 
differences, the shift q — > q — art/b enables us without any loss of generality 

1 The solution to this equation has been given by various authors with assumptions of 
even/oddness on the functions appearing flq, Ofl, u% or assumptions on the nature of B pl| ■ 
The general solution was given in ra, O]. The derivation we shall present later in fact yields 
the even/oddness assumptions of these earlier works. 

2 When b = then Vjk is just —-4-U{q), where U(q) is the potential energy of the system. 
When 6=1 and a = then the Hamiltonian for the system is pq] H = ^\- cosh 6^ Vj(q) 
where 9j is the rapidity canonically conjugate to qj and Vj(q) = Ilfe/i Vjk- hi this case (||) 
corresponds to the flows of S± =J2j e ±6j Vj(q). 



to set a = when & / and we shall do this where appropriate. Our first 
stage of generalisation then is to consider a matrix of pairwise (though in prin- 
ciple distinct) interactions and it is just this limitation to pairwise interactions 
Vjk = Vjkilj — Qk) that enables us to derive functional equations of a given type. 
Yet rather than just a single functional equation our extension allowing different 
interactions Vjk now leads to a system of functional equations, and this inter- 
play of matrix relations and functional equations appears new. Although our 
ansatz by its very form must include the Calogero-Moser and Toda models, our 
approach shows how they may be unified by the study of one functional equation. 
We find the functional equations needed to construct a Lax pair yielding the 
equations of motion (pD are of the form (for 6^0) 



<f>i(x + y) 



<p2 (x) 






4 (x) 
05 (x) 






(3) 



Elsewhere we have shown, 



Theorem 1 j^j The general analytic solution to the functional equation (|3[) is, 
up to a Q action given by (WQWi), of the form 



01 (x) 



®(x;v) ffa(x) 
$(x;//)' \M%) 



$(x; v) 
§(x\ v) 



and 



04 (x) 

05 (x) 



$(x;/i) 
$(x;/i)'7' 



Here 



*(x;,) ^^ 

a[y)a[x) 

where a(x) = ct(x\u>,uj') and C(x) = ^y are the Weierstrass sigma and zeta 
functions. 

The symmetries Q of (0) will be described in section three. The proof given in 
is constructive and the transformations needed to obtain the solutions may 
be readily implemented. 

The case b = is more problematical. In this case we obtain functional 
equations of the form 



M x + y) = M x + y)(M x ) - My)) + 



02 (x) 03 (y) 



(4) 



Certainly we may take the limit b — > to our solutions of (|3[) (which, for example, 
will give the Calogero-Moser model as the nonrelativistic limit of the relativistic 
Calogero-Moser model) to obtain solutions of (f|) but at present we don't know 



the general solution to (|]). We can show however that known nonrelativistic 
models are solutions to this equation, together with new potentials such as 

V jk (x) = aja k p'(x). 

When dj = a k this yields the usual type IV Calogero-Moser potential. 

We remark that @ and a suitably symmetrized form of (|j) are particular 
cases of the functional equation 

N 

with iV = 1 in the case 6 ^ and iV = 2 in case 6 = 0. In the case 03«+2 = 4>3i+i 
Buchstaber and Krichever have discussed (|5|) in connection with functional equa- 
tions satisfied by Baker- Akhiezer functions |14|. Dubrovin, Fokas and Santini 



18f have also investigated integrable functional equations via algebraic geometry. 

A further generalisation of the Calogero-Moser system has been to associate 

such an integrable system to the root system of an arbitrary semisimple Lie 

algebra |24|, |25[. At this stage of generalisation we have essentially the Lax pairs 



associated with A n type root systems. To incorporate more general root systems 
we may consider embedding fi : R n — > R of our n-degrees of freedom into a 
larger space with the interactions V still of the given form. We will not present 
this generalisation here. 

This paper then presents an ansatz for Lax pairs whose consistency yields 
equations of motion (^|) and the corresponding functional equations (||) and 
@. We shall show how various specialisations lead to the known systems and 
introduce some new ones. An outline of the paper is as follows. In section 
two we present the ansatz. For clarity of exposition we initially confine our 
attention to the case 5^0 returning to the 6 = case in Section five. Here 
we determine in Theorem 2 the system of functional equations to be solved, 
reducing the nontrivial equations to be solved to the form (|3|). Section three 
describes the invariances of (||) and illustrates the solution of Theorem 1 as 
a means of introducing certain elliptic function identities useful in the sequel. 
Here we will apply the general analytic solution of this theorem to our system 
of equations (Corollary 1) and then as an example show how the relativistic 
example of Bruschi and Calogero || arises. In general our ansatz leads to a 
system of functional equations and Section four looks at the constraints imposed 
on the parameters of the solution to ([]) by such a system. Theorem 3 determines 
these constraints and these are illustrated by the relativistic Calogero-Moser and 
Toda models. Further, we are able to characterise the relativistic Calogero-Moser 
model by a certain 'generic' property, Theorem 4. Section five returns to the 



6 = case deriving the appropriate functional equation. Again, several examples 
will be given. Finally we conclude with a brief discussion. 

In an earlier version of this paper we proved Theorem 1. Subsequently we 
found a direct an constructive proof which has been preseted separately || . In 
revising the present paper accordingly we have also strengthened the results of 
section 4. 

2 The Ansatz 

We shall now describe the ansatz, introducing our notation and illustrating some 
techniques useful in the reduction problem of Lax pairs to functional equations. 
Having presented the ansatz for our Lax pair we proceed to determine the restric- 
tions on the functions that appear in this. The equations we find are a natural 
generalisation of those found in ||. We will then seek the relevant functional 
equations to be solved in later sections. 

We need a few definitions in order to specify our ansatz. Let r be the fixed 
vector r — (1, . . . , 1) and denote by X d = Diag(Xi, . . . , X n ) the injection W 1 —> 
Mat(n). Further let M n = {A\A e Mat(n), A it = 0, A^ = A ij (q i - qj)} be 
a subset of matrix-valued functions of one variable. (One can also extend our 
analysis to the case of nonvanishing diagonal elements.) Note this set depends 
on a choice of coordinates {e^} with respect to which we express our matrices 
and determine the coordinate projection qi = (ej,g). A change of basis results 
in a straightforward conjugation. We denote by ej = (0, . . . , 1, . . . , 0) T the j— th 
coordinate (column) vector. Thus ej.ek = 5jk, eJ.A.ek = Ajk and X d .e^ = X^. 

Our ansatz for the Lax pair takes the formP] 

L{q) = qd + Vjar + bq) d A^J(ar + bq ) d 
M{q) = {B. [ar + bq]) d + ^(ar + bq) d C^/{ar + bq) d (6) 

where A,B,Ce M. n . When a = the matrix ^J~q d corresponds to the diagonal 
matrix D of |30]| . This Lax operator also possesses the symmetry 

L(q, b, a, A(q)) = I L(q/l, bl, a, l~ l A(l ■ q/l)), (7) 

which enables us to rescale 6 ^ to 1 and we shall later do this. It will be 
convenient to define 

G = \^ = B + -V and H = A'-C. (8) 

byj 2 



3 We remark that there is no essential change if we take 

L{q) = q d + {ar + bq) d A(ar + bq) 1 ^ 
and a similar modification to M, for any value of e. We choose the symmetrical value given. 



(This defines ^ up to a multiplicative factor that is immaterial in what follows.) 
Then we have 



Theorem 2 When b ^ the Lax pair (Qj yields the equations of motion 
and only if 



if 



with 



V, 





— Hjk + bAj k Gj k — 0, 






(9) 


Hjk — bAjkGkj + b AjkVkj — 0, 


(10) 


bAj k (Gj m — G km ) + b Aj k Vkm = b 


-rt-jm -™-mk 




(11) 


1 
1 - b 2 A jk A kj 


Aj k A k j 

A 1 A' 

A jk A kj 


= IAi n( i_ 

b 2 dx 


b 2 A jk (x)A kj 


(-*)). 


(12) 



jk 



In particular, Vj k ^ => Aj k , A k j ^ 0. 

Observe that once we have established the theorem an ansatz is entirely 
specified by the matrices A and G: given A and G then H is determined by 
@ and so C by (|j); also V, determined via A using fll2|) , together with G are 
sufficient to determine S. Thus our goal will be to find such matrices A and G. 
We remark in passing that a consequence of (0) is that V^a;) = —V k j(—x) and 
so J^ g,- is constant. We begin by establishing the theorem. 

Proof. By using the freedom to commute diagonal matrices we observe that 



L = ^(aT + bq) d (- 



b 

+ 77 



Qd 



-A+A- 



Qd 



■dA t 



H <*e £ = (^) = ([ft - ^ 
consequently yields 



or + 6g)d 2 V (ar + 6g) d (ar + 6g) d 

: [5d, A'\. The Lax equation L 



+ [q d ,A']j^(aT + bq) d . 
[L,M] 



id 



+ 



<?d 



-A + A-, 



Qd 



+ [qd,A'} 



{ar+bq) d 2y(ar+bq) d {ar+bq) d _ 

= [g d , C] + [A, (B.[ar + 6g]) d ] + A(ar + bq) d C - C(ar + bq) d A. (13) 

To solve this equation we consider the diagonal and off-diagonal terms separately. 
The fact that A, B, C have vanishing diagonal terms results in 



(Id 



(ar + bq) d 



Diag( A(ar + bq) d C - C(ar + bq) d A 



(14) 



This equation has several consequences. First, observe that if A and C have 
been determined then we obtain the sought after interaction, 



( lj 



^(a+bq j )(a+bq k )V jk (q j -q k ) = (a+bqj)\A(aT+bq) d C-C(aT+bq) d A 



k^j 



j j 



That is 



A k Au- 
Vjkilj — Qk) — Aj k Ckj — CjkA k j = n 3 „ J = —Vkj(qk — Qj)- (15) 

^jk ^kj 



The second and crucial point is that (|T4]) reduces the off-diagonal terms in our 



consistency equation (|13|) to a linear equation in ar and bq. Utilising our defini- 
tions of G and H we find the off-diagonal terms of ( |i3|) yield (for j ^ k) 



[q d ,H} + [(G.[aT+bq]) d ,A\+bA{V.[aT+bq]) d +C{ar+bq) d A-A{aT+bq) d C) = 0. 

/ jk 

When b ^ this equation reduces to considering 

([q d , H] + 6[(G.g) d , A] + b 2 A(V.q) d + &Cg d A - bAq d c) , & = (16) 

as the r terms of the previous equation are reproduced by taking q oc r. We will 
consider separately the 6 = situation which gives the equation 

([q d , H] + a[(B.r) d , A] + a[C, A]) jfe = 0. (17) 

We may solve the linear equation ( JT5| ) by taking particular choices for q. 
Substituting q = e k , q = tj and q = e m {m ^ j, k) in (pT6|) yields the equations 
d), (0) and (0) respectively 

We have thus shown how the Lax pair determined by the data A, B and C 
fixes the potential V via (15) and consequently matrices G and H such that (§), 
(|T0| ) and (|11]) must hold. It remains to show that V is also given by (|I2|). First, 
if Ajk = then by @ and (H) we find Cjk = and so by (^) V^ = as well; 
thus ( ]T^ ) holds in this case. Finally, if Aj k ^ 0, then upon adding ([T^) and (Q) 
we obtain 

bV kj (x) = G kj (x)-G jk (-x). (18) 

Making use of the antisymmetry of V together with the definition of G gives as 
a consequence of (|18|) that Bkj(x) = Bjk(-x). Equation ([12]) now follows from 
(|15|) after making use of this symmetry and the expression 

C jk = A' lk - bA jk G ]k} (19) 

which follows from (^j. Therefore we have established (|L?]). 

The converse of the theorem follows from our initial remarks that the matrices 
A and G together with the definition of V suffice to determine the Lax pair. 

□ 



Our task therfore is to find matrices A and G for which (|9|-|T2|) are satisfied. 
In order to understand this system of equations it is helpful to consider the 
consequences of an entry of A either vanishing identically, or otherwise. 

Lemma 1 If Ajk = for some j,k then Hj k = Cj k = Vj k = 0. Further, for 
any m ^ j,k for which Cj m , C mk ^ there exists a constant aj m k such that 
A jm (x) = a jmk C jm (x) and A mk (y) = a jmk C mk (y). If a jmk ^ then 

A jm (x) = a jm e x/a J mk tfjjmix), A mk (y) = a mk e y,a ^ k tp mk (y). (20) 

for constants ctj m , a mk . 

Proof. The vanishing of Hj k , Cj k and Vj k is immediate. Further our assumption 
means that the right-hand side of (|ll]) vanishes. Because Cj m and C mk do not 
vanish identically, flT5| ) entails that neither do Aj m and A mk . Thus the first row 
of the determinant must be proportional to the second and we have the second 
assertion of the lemma. Upon making use of flT§|) we obtain 



Aj m {x) — aj mk {A- m {x) — uAjmlxjGjmyx)) — aj mk {A- m {x) — Aj m {x) 



^jmO)- 



Yjm \X) 

When aj mk ^ this gives 

Aj m {x) _ 1 _ j^mW 

-™-jm\X) Qjmk YJm\X) 

and the first part of (p0|) follows upon integration; the expression for A mk (y) 
follows similarly. 

□ 

Lemma 2 If Aj k ^ for some j,k then V k j(x) = G k j(x) — Gj k (—x) and 
B k j(x) = Bjk(—x). Further, for any m ^ j, k: 

1. if A mk = 0, then Gj m = G km = c 1; a constant; 

2. if A mk ^ and Aj m = then Gj m = G mk = c 2; a constant; 



3. if A mk ^ and Gj m — G mk ^ (and consequently Aj m ^ 0) then pi| j may 
be written as 



Aj k 



A ■ A , 
^jm - rl mk 

^ jm ^mk 



Hj rn /\ rnk "T 



A ■ A , 

■"-jm - rl mk 

A' A' 

.im rnk 



jm ^mk *-* jm ^mk 



(21) 



Proof. We have already proven the first part of this lemma in our discussion of 
the theorem, ([Tj|). Now if A mk = for some m, then V km = by ([12]) and so 
only the first term of ( |Tl"D is nonvanishing. This means Gj m = G km and as these 
are functions of different arguments, they must be constant. 

Assume A mk ^ 0. Then employing fll8|) for these indices enables us to rewrite 



b{Aj k — b Aj m A m k){Gj m 
The remaining cases are now straightforward. 



A ■ 
A' 

jm 



A 



rnk 



A' 



rnk 



(22) 



D 



Clearly (|21|) yields a nontrivial relation between the desired functions A and 
G and we shall now show that this equation may be recast in the form ([|). With 
this aim in mind it is convenient to rewrite fl2~T|) in the form 



A jk {x + y) = A jm (x)A mk (y) - ip jm (x)ip mk (y) 



"■jm\~ 



A m k(y) 

A'jm( x ) A' mk (y) 

1pjm{x) 1pmk{y) 
Vjmtx) fmkiv) 



(23) 



where A = b A. As we have already remarked, when & ^ we may rescale so 
that b = 1 and for the remainder of this section and until section 5 we assume 
that we have done this. 

Lemma 3 Equation $2l\ ) may be rewritten in the form (fjj where either 



A jk (x + y) 



Aj m [X) [1pj m ( X 



Ajrn(y)/lpjm(y) 

A mk (x) / j) mk (x) A mk (y)/ip mk (y) 
l/ip jm (x) l/ijj jm (y) 

l/4>mk(x) ±/4>mk(y) 



(24) 



or 



A jk (x + y) 



c 2 A mk (x) c 2 A mk (y) 

^mk{ x ) ^mk{y) 



G mk (x) 

1 



G m k{y) 

1 



for some constant c 2 according to whether 
not. 



jm 



. X 



1pmk{x) 1pmk{y) 



^jm{y) 



(25) 



7^ vanishes or 



Proof. Adopting the shorthand d = 

C -^-jm-rt-mk 



dj + d k observe that 



.4 



A 



A' A' 



■jm 



ink 



]» 



rnk 



Thus we may rewrite fl23|) as 



-^■jk Sijm-^-mk tyjmtymk 



O Aqrn,A 



jm-rt-mk 



(26) 



Now because Ajk - 
sides of (EBl) shows 



dlpjmlpmk 

Ajk(xj — Xk) we have dAjk = 0. Thus applying d to both 

" -^-jm-^mk 



d 



0. 



-t, i p Jm mk — A ■ At- 



d^j m 1pmk 

and consequently the ratio here is a function of x 

Xk)- If we set x = Xj — x m and y = x m — Xk then (|26|) takes the form 

A jk {x + y) = A jm (x)A mk {y) - ^jm{x)i) m k{y)A jmk {x + y). 



(27) 



The left hand side of this equation is symmetric under the interchange of x and 
y. Performing this interchange and subtracting the resulting equation from fl2T| ) 
shows 







■fi-jm \X) 
■fi-mkyE) 



Ajmiy) 
A m k(y) 



I'm 



. X 



^mk{x) 



^jm(y) 
i>mk(y) 



A jmk (x + y). 



(28) 



Two possibilities now arise, each leading to an equation of the form (Q). 



Suppose first that 



1pjm(x) 1p jm (y) 
^mk(x) 1pmk(y) 



A jmk {x + y) 



^ 0. Then 



-rt-jm \X ) 

-™-mk\X) 



Ajm{y) 

Amkiy) 



tpj m {x) 

^Pmk{x) 



^mkiy) 



Upon substituting this expression into (|27| ) one obtains 

1pjm(x) 1pjm(y) 
^Pmk(x) ^mk{y) 



In the case when 



some constant c\. Therefore G 



jm 



X) 



■™-jm\X) 



C2A mk (x) and hence Cj m (x) 



- we have ipj 
G m k(x). Likewise from 

— tyGmkiX 



X) 



Ci^mk(x) for 
we have that 



into ( pl|) obtains fl2~5|), again of the stated form. 



Substituting these relations 

□ 



3 The Functional Equation 

Thus far we have reduced the consistency requirements for the Lax pair (||) to a 
functional equation of the form (Q) and this section looks briefly at this equation. 
Particular cases of this equation have been described in the literature (7], |6| . 



10 



Here we describe the invariances Q of (|3|). Theorem 1 gives a representative 
of each Q orbit on the solutions of (|3|) with a particularly nice form. For later 
calculations it is instructive to see how the stated solution satisfies (^), an exercise 
involving some elliptic function identities. We end the section by deriving the 
relativistic Calogero-Moser model found by Bruschi and Calogero |J. 

First observe that a large group of symmetries Q act on the solutions of @ . 
The transformation 



(f>i(x), 



h{x / 
clearly preserves 






e Xx (p2(x] 



^^m^j-iv 



e A x 4>4;(x] 
e x " x Mx) 



provided 
A + A' + A" = 0, U,VeGL 2 , 

Further, (H) is also preserved by 



and 



detU = c det V. 



and 



<px(x), 



h [x 



M x , 

03 ( x ) 



M x )\ (M x ) 

03^)/ ' \M X ) 



04 (x) 
(f) 5 (x) 



4>±(x] 
fa (x) : \M x l 



M x ] 

fo( x ] 



^<'<SV<!:; 



(29) 
(30) 

(31) 
(32) 



These symmetries enable one to find a solution of (|J) on each Q orbit with 
a particularly nice form. It is instructive to see how the stated solution satisfies 
(131). ^From the definition of the zeta function we have 



Thus 
$(x;z/) <&{y;v) 



( ln$(z; v))'= -t(v -x)- C(x) + Ciy). 

= ®(x;v)®(y;v) (ln$(y; v))'—(ln$(x; v))' 

= $(x; v)${y; v) V{y - x) + ((x) + ((-y) + ((y - v) 



Upon using the definition of $ the right hand side of this equation takes the 
form 



$(x + y;v] 



a{y — x)o{y — y)cr(x + y) 



<j(v - x - y)a(u)a(x)cr(y) 
After noting the two identities |34| 



C(v-x) + ({x) + ((-y) + (( y -v) • (33) 



C(x) + C(y) + C(z) - ({x + y + z) 



a(x + y)cr(y + z)a(z + x) 
a(x)a(y)a(z)a(x + y + z) 



(34) 



11 



and 



p{x) - p(y) 



a(y - x)a(y + x) 



(35) 



a 2 (y)cr 2 (x) 

we find ( |3~3"D simplifies to $(x + y\ v)[p(x) — p(y)], where p(x) = —('(x) is the 
Weierstrass jp-function. Putting these together yields the addition formula 

$(x;v) $(y;v) 
$(x;v)' $(y;i/)' 



P{x) - p{y) 

and consequently a solution of @ with the stated form. 
A consequence of Theorem 1 then is 



(36) 



Corollary 1 If Aj k satisfies ft2l\ ) then 

Ajk{x) = Cjk 



$fo Vjh) c x ikX 



^ K $(x,/i iyk 
for some constants Cjk, Vj k , f^jk an d Xj k - 



(37) 



Example The relativistic example of Calogero and Bruschi arises as a particular 
case of our ansatz when a = 0, 6 = 1 and 



A jk (x) = (1 - 5 jk )a(x), B jk (x) = (1 - 5 jk )(3(x), 
In this case (E3) takes the form 



C jk (x) = (l-6 jk )j{x). 



a(x + y) 



a(x) a{y) 

T](x) r](y) 

1 1 



(3* 



where r)(x) = (3(x) + \v(x). Comparison of 
shows the solution to be given by 

&(x:v) 



a(x) 



7](x) 



ln$(x;yu))'= - 



with the general solution of 
1 pf(x) - p/(fi) 



$(x;/i) '" 2 p{x) -p(fi) ' 

In this case a(x)a(— x) = (p(x) — p(u))/(p(x) — p(fj)) and so by ( |T2"D 

p/(x) 



V[X) 



p(ix) -p{xY 
and we have recovered! 3 ! the results of references |E9| and 



4 We remark in passing that the the scaling of the elliptic function 

p(tx\ttu,tuj') — t~ p(x\ui,uj') 

is accounted for by the scaling of the Lax operator (pi) and doesnt appear to give any new 
potentials. 
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4 Application of the Functional Equation 

Thus far we have discussed the functional equation (^) in isolation while our 
application involves a system of such equations. In solving this system we en- 
counter further constraints. To see how these arise consider ([24]). Our theorem 
says 

A >Ax) = Cju ) ' 3hJ e x '" x 
$(x;/i jfc ) 

(for some constants Cj k , Vjk, Pjk and Xj k ) and further that 

A jm (x)/i) jm (x)\ = , , _ Wxu ( $0, v jk ) 

A mk {x)/ip m k{x)J J{ ' \&(x,u jk ) 

V^jm{x) \ = f(x)e \"xy ( Qfafijk) 
l/^mk{x)J n ' \<$>'(x,H jk ) 

for an appropriate function f(x) and matrices U, V such that 
X jk + X' + X" = det U = c jk det V. 



Now if Aj m and A mk are also given by fl37|) we encounter restrictions on the 
possible parameters appearing: 



Theorem 3 Let Aj k , Aj m and A mk have the form $3% ) and be related by ([£ 



Then the constants determining these solutions are related by 

Vjk r'jk "jm pjm ^mk Pmk v / 

\k + C( U jk) - C(Pjk) = \m + C(Vjm) - C(Pjm) = X mk + ((v mk ) - ((Pmk) (40) 

and c i k = ^T) T i k > c i m = ^S~) T i m ' Cmk = $^) Tmk where the T ' s sat%s fy 

Tjk _ Cr[yj m + l> m k ~ Vjk) 



'jm'mk *~>\Pjm ~T~ Pmk pjk) 

Letting fij k — Vj k = c and Xj k + C{ u jk) — ({Pjk) = P then 



(41) 



cr(c + Vjk ~ x) 



O-jVjm ~ X) ^ px _ O-(Vmk-x) ,, 

jm ' jm / \ " i s\mk 'mk / . \ 

a(c + v jm - x) a{c + v mk - x) 



and similarly 

■"■jm 1~j m - f ; ~e , -rt-mk T~mk 7 I \^ 

Finally 

Gjm(x) = ~((X - fXjm) + C{X + Hmk ~ Pjk) + COUSt (43) 

and similarly for G mk (y) with the same constant appearing. 
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The proof of this theorem is rather lengthy, making repeated use of the ellip- 
tic function identies introduced in the previous section; it is given in appendix A. 
The first two relations ( |55| ) and ( JfOp follow by equating poles and zeros amongst 
the various terms while the relation (flip comes from the determinantal con- 
straint. The final constraint fl43|) arises by considering (p2|) , which may be recast 

as 

Aj m {x)A mk {y) 



Gim(x) — Gmkiv) — dial 1 — 

We remark that when A jm (x) = c 2 A mk (x) (and so u jm = u mk , fi jm = fi mk 
Ajm = A m fc, Cjm = C2C mk ) several of these relations are immediately satisfied. 

It is worth reflecting on this theorem. Given any three Aj k , Aj m and A mi 
of the form ( |3?| ) and connected via (2T), the constants determining these func- 



tions are restricted. In particular, suppose every entry of the Lax matrix A is 
nonvanishing and Gj m — G mk ^ for every triple of distinct indices. Then the 
theorem holds for every triple Aj k , Aj m and A mk . Consideration of fl43|) shows 



that if this is to define a function Gj m then 

/J>jm + Mmfe - l^jk = H (44) 

for some fixed constant \x and every distinct triple j, m, k. Now (fP|) holds for 
every distinct triple if and only if (for each j, k) 

fJ'jk = A* + 1^3 ~ A*fe 

and similarly Uj k = v + Vj — v k . In this case Tj k = a(fj,)/a(u) and so 



A jk (x) 
That is 



o{li)a{v + v 3 -v k -x) ^ px 
<y{v) a(fi + Vj -v k -x) 



A- k (x-- x k ) = ® {Xj - Vl - {Xk - Vk)] u K lp - CM+CM]x 
1 J $(xj-Uj - (x k -v k y,n) 

and we have recovered the relativistic Calogero-Moser model described in the 
last section. Now Gj m = G mk <=>■ Gj m = G mk = constant. Thus we obtain the 
following 'generic' description of the relativistic Calogero-Moser model: 

Theorem 4 A Lax Pair of the form (^) for which the matrix A has no van- 
ishing entries and for which the matrix G has nonconstant entries describes the 
relativistic Calogero-Moser model. 



We have just considered the situation where every entry of A satisfying 
We conclude the section by considering the opposite extreme where no entry 
does. 
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Example Here we adopt the ansatz 



A 



/0 oi ...0 \ 
1 a 2 

1 1 o 3 00 



1111 
\1 1 1 1 



O n -l 

1 / 



and where we shall determine the a^s 7^ 0. We will work through the various 
cases determined by Lemma 1. 



1. j < k — 1. Then Ajk = and the only possible nonzero term in (|i~TD is 





Au+i A j+lj+2 



W7+I ^ 



j+lj+2 



Thus for each j < n — 1 



C,,' 



ii+i 



C.7+1 



i+ii+2 



A/j+i A j+lj+2 
and consequently for each j < n 



v 

ft- ; 



ii+i 



.4 



G 



ii+i 



ii+i 



A. 



(45) 



2. j = k — 1. In this case A/^+i 7^ and now if 

(a) m < j then A mfc = and G jm = G j+lm . 

(b) j + 2 < m then G jm = G j+ i m . 

(c) m = j + 2 then G^+2 — Gj+ij+2 + Vj+ij+2 = whence upon using 



(0) 



Gjj+2 - G j+ i 



& 



i+ii+2 



•j+ij+2 



1-A 



3. k < j. Then A^ = 1. Now fl22"|) becomes 

and we find 
(a) if m < k — 1 or k < j < m — 1 that Gj m 



j+lj+2 



A ■ A , 

/I' /I' 

jm mk 



^mki 



(46) 
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(b) m — k — 1 



(c) j = m — 1 



Gjk-i — G 



k-lk 



Gjj+i - G j+1 



j+ik 



A' 
^k-lk 

1 — Ak-ik 



A' 
1 ~~ A jj+i 



(47) 



(48) 



with no constraints arising when k < m < j. Further, when k < j + 1, we 
have Vjk = = Gjk — Gkj and so 



Gjk = G kj \j - k\ > 1. 



(49) 



Now case (2a) tells us each column of the matrix G is constant below the 
diagonal while (2b) tells us each column above the superdiagonal is constant. 
Combining this information with (H9j) enables us to parameterize G as 



G 



/0 9i 





. 


.. 


^ 




(0 


1 


1 1 . 


• l l \ 


1 


92 













1 





1 1 


1 1 


1 1 





9s 








+ d 


1 


1 


1 


1 1 


1 1 


1 


1 





9n-l 




1 


1 


1 1 


1 


\1 1 


1 


1 


1 


J 




u 


1 


1 1 


1 oj 



where d is a constant. Comparison with ([45]) -([48]) shows we are left with two 
equations, 



9j 



= V, 



ii+1 



-1 - g. = A + d, 



with solution 



aj(x) = 1 



1 + c jj+ ie( x+d > 

and 9j(x) = — (A + d)dj(x). Here Cjj+\ is a constant that may be removed by 
shifting the x's. 



We have recovered the relativistic Toda lattice of [ 30(1 - O ur construction has 
given the Lax pair (b — 1, a — 0) L — yfqd (I + A)^fqd and 



i9iQs 







0' 



M = {B.q) d +^~ d C* 



-dL+{d^q)I+\ ° \gi<ii + \92q-i 

3=1 V 
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+- 






-9i 




. 

~92 


.. 
























~9n 





/ 

Upon denning a conjugate Lax pair L N , M N by L = \fq~dLN 1/ ' \fq~d and Mm 
^/Vq~d(M + dL — (dJ2 r j=i<l)I)\/<ld+ VQd/VQd we obtain the Lax pair 



L 



N 



M 



N 



When A = 1 and d = this is the Lax pair of 



5 The Nonrelativistic case. 



fqi 
Qi 
Ox 


a x q 2 
<l2 a 2 q 3 

q.2 q-i 











q\ 

\qi 


q.2 q-i 
q.2 q-i 






0"n-lQn 
q n J 


( 9ih 




~9ih 

92<l3 - 




-92C(3 




° \ 






V o 












-9n-iq n 
/ 



Our discussion has thus far focussed on the b ^ case of our ansatz for the Lax- 
pair and this has yielded the relativistic Calogero-Moser and Toda systems. We 
shall now consider the case 6 = and see that this includes the nonrelativistic 
limits of these systems. The nonrelativistic dynamics is less constrained than 
the relativistic situation and this reflects itself in a more complicated functional 
equation which we have not been able to solve in general. In this section we 
shall first obtain the relevant equation and show how it encompasses the non- 
relativistic Calogero-Moser and Toda systems as well as that of Buchstaber and 
Perelomov 
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Our approach leads to a new derivation of the Calogero-Moser 
model in which we also determine the various symmetry properties of the func- 
tions entering the ansatz. Recalling that the diagonal entries oi A, B G M. n 
vanish we begin with: 

Theorem 5 Equation flll[ ) yields the functional equation 

J2 A jk (B jt - B kl ) + A! jt A lk - A fl A' lk = 0. (50) 
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Further, each of the functions *&ji k = Ajk(Bji — B k {) + A'^Aik — AjiA' lk appearing 
as the terms of this sum, depend only on the combination qj — q k . Thus with 
x = qj — qi and y = qi — q k we have 

Vjikix + y) = A jk {x + y){B jl (x)-B k i{-y))+Af jl (x)A ik {y)-A jl {x)A! lk {y). (51) 

Proof. In order for ( fT7|) to remain true for all q we must have 

H jk = (52) 

and therefore Cj k = A'- k . In fact this is the h — > limit of equations @ and 
(|10|). With these simplifications ( |17D then becomes 



Ajk{Bj k — B k j) + 2_^ A jk (Bji — B kl ) + A'^Aik — A^A\ k = 0, 



i 



and with the conventions stated above we have 

Now the only dependence on qi (I ^ j, k) in fl5"0|) comes from the Z-th term of 
the sum. Thus upon taking the derivative di of this equation we see that 



= diiA^Bfl-B^+^Aik-AflA!^ = -{d.+d^A^B^-B^+A^-A^ 



ik j 



from which we may conclude that ^ji k = ^jik{<lj — Qk) as stated. We remark 
that the b — > limit of equation (flT|) divided by b gives the quantity \& j mk . □ 



Corollary 2 Solutions of fl5l\ ) satisfy the functional equation (||) with N = 2. 
Proof. Upon interchanging x and y in (|5l"D and subtracting we obtain 

A jk (x + y) 



'ji(x) +B k i(-x) 
1 


Bji{y) + B kl (-y) 

1 


+ 


A'jtix) A'^y) 
Ai k {x) A ik (y) 


(53) 


+ 


A' lk (x) A' lk {y) 
Aji(x) Aji{y) 


= 


(54) 












(55) 



Simple rearrangement of this gives equations of the form (^) or of (|5|) with 
N = 2. 

□ 

Before further investigating (^Dj) it is instructive to see how the nonrelativis- 
tic Calogero-Moser and Toda systems arise in this context and present a new 
example. 



Example The Calogero-Moser reduction follows when we assume the functions 
Aj k and Bj k dont depend on the indices j,k. Upon setting Aj k (qj — q k ) = 
A(qj — q k ) and Bj k (qj — q k ) = B(qj — q k ) the function tyj lk defined above takes 
the form 

tyjOs = A(q j -q k )(B(q j -q l )-B(q k -qi))+A'(q j -q l )A(q l -q k )-A(q j -q l )A'(qi-q k ). 

Now the lemma asserts that ^j lk is independent of qi (I ^ j, k) and so (for 
example by setting qi — — qi>) we see that ^ji k = ^ji> k for I, V ^ j,k. Thus 
( |50| ) takes the form 

= A(x + y)(B(x + y) - B(-x-y)) + 

(n - 2)(A(x + y)[B(x) - B{-y)\ + A'(x)A(y) - A(x)A'(y)), (56) 

where x = qj and y = —q k . Letting B denote the odd part of B, then upon 
interchanging x and y in fl56D and adding leads to 



(n — 2) 
= A(x + y)(B (x + y)+ K -^-[B {x) + B (y)}). 

Now the functional equation (n > 1) 

(n — 2) 
= B (x + y)+ K -^[B (x) + B (y)) 

only has B = as a solution and so if A(x) ^0 we may assume 5 in ( pB] ) is an 
even function. In this case the leading term vanishes and we are left with 

= A(x + y)[B(x) - B(y)} + A'(x)A(y) - A(x)A'(y) 



which is the functional equation (jl|) obtained by Calogero [|16|] and that has been 
variously solved 0, ^6[ |2^, 0, |IT| . Observe that by setting y = —x in this 
equation and using the fact that B(x) is an even function we may deduce that 
A(x) is an odd function and so we have obtained the symmetries of A and B 
normally (16], |26| imposed as constraints on the Lax-pair. 

Example The Toda reduction follows by assuming Bj k and ^ jik vanish for all 
possible distinct indices. From this we deduce that 

A jk (x) = a jk e x J kX 

for some constants aj k and \j k . The vanishing of ^/ji k relates these constants 
(for j ^ l=£ k) by 

(Aj/ - Mk)®jia lk = (57) 
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while the nonvanishing of an interaction Vjk means (using (|I5|)) that 

(Ajjb - \kj)®jk®kj ^ 0. (58) 

We thus seek solutions of ([H]) and fl5"5|). Observe that if Vjfj and Vj*. are both 
nonvanishing then Xji = \i k ^ \ ki = Xij. Consequently if Vlj, Vji and V\k are 
nonvanishing then (for example) Vjk = 0. To see this suppose to the contrary 
and note that ^7. = means \j = Xjk, while ^ jki = means Xjk = Xm and so 
Xij = Xm- However our previous observation shows Xij — . . . — Xik 7^ Xm, yielding 
a contradiction. Therefore Vjk = 0. Our ansatz means we can only achieve 
(possibly cyclic) chains of nonvanishing interactions, V\ 2 = V23 = . . . = V^_i„ 
(= V n i). Setting a ii+ i = 1 = a i+ u and X ii+ i = —X i+ n say, with all other 
oiij and X^ vanishing, yields the (periodic) Toda lattice. Equally the work of 
|20f shows how to obtain the Toda systems as limits of Calogero-Moser models. 
Example By taking 

Ajk(x) = $(2;; v) a k , Bjk(x) = p(x) a k (59) 

where the a& are constants, we find 

®(x;v) $(y;v) 
$(x;z/)' $(y;u)' 



^jik(x + y) = $(x + y; v) (p(x) - p(y)) a k a t 



a k ai = 0. 



Here we have used ([36]) to show the vanishing of ^ jik- By Theorem 5 we therfore 
have a new Lax pair associated with the potentials (|2]) 

Vjk(x) = ajakp'(x). 

When cij = ak this yields the usual type IV Calogero-Moser potential. 



Although we cannot as yet solve (|50"D or ([51]) in general, we can say a little 



more according to whether AjkiBji — Bki) = or Ajk 7^ 0. For the first we note: 

Lemma 4 The functional equation 

F(x + y) = <t>{x)i)'{y) - <j>' {x)ij>{y) (60) 

has, up to symmetries, the solution F(x) = ciCsexp(Xx), 0(x) = ciexp(Ax), 
ip(x) = (02 + csx) exp(Ax). 

This equation is a particular case of the equation 

f(x + y) + g(x-y) = Y^ fj(x)9j(y) 

which has a long history ||; the solution, along standard lines, is given in Ap- 
pendix B. In our context it yields 
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Corollary 3 If Aj k (Bji — B k i) = the solution to ffi]\ ) takes the form ^ ji k {x) = 
(ciC 3 -c 2 c 4 )exp(Aa;) ; A fl {x) = (c 2 +c 3 x) exp(Ax) and A lk (y) = (c 1 +c i y) exp(Xy), 
where c 3 c± = 0. 

Remark: The Toda reduction was a particular example of this corollary. The 
various constants appearing are not all independent, but related by (p0|). 

We now suppose that Aj k 7^ 0. Upon setting 

A' A , — A ■ A' 
Jfc_2 - A~ k 

for j ^ k and $.,%,• = we may rewrite flSDp as 

^ (5 rm - B fem ) + % fe = 0. (61) 

Lemma 5 For each j , k and I for which Ajk, Am and Aij are nonvanishing then 
( ]3Z| j yields the functional equation 

® jk + $ w + $ y = (62) 

and consequently 

$ ifc + $ fci = (63) 

Proof. This follows from flBTD as 

D 

Example In the case when n = 3, equation (|62[) reduces to the functional 
equation of Buchstaber and Perelomov jl3| , 



( /(x) + g(y) + &(*) ) 2 = F{x) + G(y) + H(z) x + y + z = 0. (64) 

The equation is related to the factorization of a three-body quantum mechanical 



ground-state wavefunction |31], [32|, |I7|]. In this case there is a unique m 7^ j,k 

such that 

A' A , — A ■ A' 

-j- _ J ^jm^ l mk STjm^mk 

jk ~ A~, ' 

Set x = q 2 — q 3 , y = q 3 — q±, z = q 1 — q 2 (and so x + y + z = 0) and 

f\x) = -A 23 {x)A 32 {-x), g{y) = -A 31 (y)A 13 {-y), h(z) = -A 12 (z)A 21 (-z). 
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Then 



$21 + $13 + $ 



-32 







leads to 

A 23 A 31 {A 13 A' 32 -A[ 3 A 32 )+A 31 A 12 {A 21 A' 13 -A' 21 A 13 )+A 12 A 23 (A 32 A' 21 -A' 32 A 21 )=0 

and consequently 

g(y)A 23 A 32 -f(x)A 31 A[ 3 +h(z)A 31 A[ 3 -g(y)A 12 A' 21 +f(x)A 12 A 21 -h(z)A 23 A' 32 = 0. 

(65) 
Similarly from $12 + $23 + $31 = we obtain 

f{x)A l3 A , 31 -g{y)A 32 A' 23 +h{z)A 32 A 23 -f{x)A 2l A' l2 +g{y)A 21 A 12 -h{z)A l3 A 31 =Q. 

(66) 

Now 

/'(*) = A 23 A' 32 -A 32 A' 23 , g' (y) = A 31 A' 13 - A 13 A' 31 , ti(z) = A 12 A' 21 -A 21 A[ 2 , 
so upon adding (^) and fl6"6|) we obtain 

g(y)f(x) - f(x)g'(y) + h(z)g'(y) - g(y)h'(z) + f(x)ti(z) - h(z)f(x) = 0. (67) 
Equation (1671) may be rewritten as 



1 1 1 

f(x) g(y) h(z) 
/'(*) g\y) ti{z) 



0. 



x + y + z = 



(68) 



which is the equation studied by Buchstaber and Perelomov. The solutions to 
(|5"8D are called nondegenerate if each of f(x), g(y) and h(z) have poles lying 
in some finite domain of the complex plane. Degenerate solutions may then be 
obtained from these. The nondegenerate solutions to (|68| ) are given by 



f(x) = ap(x-ai)+/3, g(y) = ap(y-a 2 )+ (3, h(z) = ap(z-a 3 ) + /3, (69) 

with a\ + a 2 + 03 = 0. 

We now observe that, although we have not yet specified A\ 2l A 23 and A 3 \, 
we have in fact obtained the interactions in the present situation. We have 



and so 



Vjk — —(Aj k A k j) 
V jk{Qj - Qk) = -®p'{qj ~qk + aj - a k ). 



(70) 
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Thus our functional equation determines the interaction for us. As for the Lax- 
pair we may simply choose Aj k (x) = —A k j{—x) or some other form that suits our 
purpose. Using the addition properties of the elliptic functions another choice 
for Aj k (x) could be 

A jk {x) = Vtt^TTT \ ,\ \ 71 

a(b)a(x - Xj + X k ) 

where ap(b) = —(3 and a« = Xj — X k for cyclic i,j, k. 



6 Discussion 

This paper has introduced a new ansatz (|) for a Lax pair describing systems of 
particles on the line interacting via pairwise forces (0). Unlike existing ansatz 
we allow these forces to depend in principle on the particle pair, and so the one 
ansatz encompasses for example those of the Calogero-Moser and Toda systems 
within a unified framework. A consequence of allowing varying pairwise inter- 
actions is that the consistency equations for the Lax pair now become a system 
of functional equations rather than a single functional equation. Our approach 
has been to first study the constituent functional equations, of interest in their 
own right, and then to examine the contraints imposed by the system of which 
they are a part. 

Two quite interesting functional equations (|3|), (H) arise in this manner. The 
first, which arises when b ^ 0, has a large group of symmetries acting on it and 
we have been able to give its general analytic solution with appropriate orbits 
corresponding to the relativistic Calogero-Moser and Toda interactions. It is 
this large symmetry group of (H) that enables us to relate previously distinct 
functional equations and different physical models. We remark that a particular 
case of this equation has recently arisen in the work [12| (see their equation (13) 
and Lemma 10) which examines the connection between functional equations 
and Dunkl operators. 

Unfortunately we have not been able to say as much about the functional 
equation (f|) or the associated system of equations when 6 = 0. Certainly the 
b — > limit of our general solution yields a b = solution, corresponding to an 
appropriate nonrelativistic limit, but the nonrelativistic equations are less rigid. 
Similarly we note that both (|J), (|55|) as well as the functional equations satisfied 
by Baker- Akhiezer functions |TJ| are particular cases of (El). 

The final step in our approach has been to examine the contraints imposed by 
the system of functional equations on the parameters appearing in the solutions 
to ([3D and (H). The constraints for the relativistic system were encapsulated 
in Theorem 3. Although we have a conceptually straightforward unification of 
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various ansatz for Lax pairs, this stage of our approach is the most tedious as it 



can often involve case by case analysis. We plan to return to the equation (fy]) in 
a later work. Finally we have shown how the examples of the known relativistic 
and nonrelativistic Toda and Calogero-Moser models arise in our approach as 
well as introducing a new system. 

Acknowledgements: One of us (V.M.B.) thanks the Royal Society for a 
Kapitza Fellowship in 1993 and the EPSRC for a Visiting Fellowship. 

A Proof of Theorem 2 

In proving this theorem we shall consider the constraints imposed by (|2"1|) and 
(|25|) separately. In both cases we proceed by first finding relations amongst the 
constants appearing in ( p7|) by equating the poles and zeros of the various terms 
given by our first theorem. This determines the functions ipj m and ip mk up to 



the action of the function f(x) and an exponential. By comparing with (|22|) we 
then determine Gj m and G mk up to constants. 



A.l Consistency for ([ 

We begin with (JM[). In this case the theorem says 

A jm (x)/ip jm (x) = f(x)e~ x ' x (u u ®(x; u jk ) + u 12 $'(x; u jk )) 

and 

lMmO) = f(x)e x " x (v u $(x;fi jk ) + v 12 &(x; fi jk )). 

The ratio of these two equations when 



-™-jm\X ) Cjr, 



&(x: V; 



' u 3 m > „A,- m a; 



Sf^X, flj r 



']m J 



yields 



$(z; u jm ) , x _ mi$(j; Uj k ) + u 12 &(x; u jk ) , x 
$(x; (j, jm ) vu$(x; n jk ) + v 12 ®'{x; fx jk ) 

= u 12 $(s; u jk ) [u u /u 12 + C(^-fc) - C(x) - C(^-fc - x)] cX . kX 

v 12 $(x; /i jfe ) [i)n/i;i2 + C(/^fc) - C(^) - C(/^fc - »)] 

(We will assume u± 2 and f 12 are nonvanishing and later see that this is this case. 
Certainly by considering the behaviour of ( [72|) as x — » we see iti 2 = <^ f i2 = 
-v=> A, m (x) = cA-fc(i) e A;c .) Now the left hand side has a zero at Vj m and pole 
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at fij m . Equating these with the right hand side (for Uj m ^ Vj k and \ij m ^ /ijk) 
shows 

«ll/«12 = ({Vjk ~ Vjm) ~ ({Vjk) + ({Vjm) 
Vu/Vja = CiPjk ~ fijm) ~ ((Hjk) + CifJ'jrn)- 

Thus, after making use of (|34"1) , 



u u /u 12 + ((u jk ) - C(x) - C{U jk - x) = J J — J (73) 



and 



vu/v 12 + CC/j-ja) - C(») - C(a*j* - ^) 



v(vjm)(r(vjk - v jm )a(x - u jk )a(x) 
a(fi jk )a(x - fi jm )a(x - fi jk + /ij 



Utilising ([TUP and ( |7i|) in ( |T2| ) now gives 



• (74) 



(Aim-KtemJ-Cfoym))* _ M 12 ofo - l/jfc + Vjm)<T{(J, jk ~ fijm) J\ nk+ C(^k)-C(^k.))x 



^m 



u 12 er(a; - /ijfc + Hj m )a(p jk - u jm ) 



from which we deduce 



^jm + C{ U jm) — C(^jm) = ^jk + C{ U jk) ~ Ci^jk) 
CjmV 12 = U 12 . 



^From ( [76|) and our expressions for m and — we find 

"jm ^jk 



W2 C12 

^11 Mil 
l>12 Ul2 



(75) 
(76) 
(77) 



(78 



The same considerations now applied to v4 m fc(x)/?/> m fc(a;) rather than Aj m (x)/ipj m (x) 
similarly show 

U2l/U 2 2 = C{Vjk - Vmk) - C( U jk) + Ciymk) 

^2l/^22 = CiUjk ~ Hmk) ~ CiUjk) + C(^mfe) 

Vjk — v mk = Hjk — l^rak (79) 

A m fc + C( z/ mfc) — C(/ i mfc) = Ajfc + C( z/ jfc) — Ci^jk) (80) 

C m fc^22 = «22 (81) 

^mk ~ ^jk = • (82) 

V 2 2 «22 



Combining these relations with (|75|) and ([76D thus proves (J39|) and ([40]) for the 
case being examined. 
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We have yet to impose the constraint det U = Cj k det V. Now 

det U _ U12U22 (^t _ W 



Cjk 



detV vi 2 U22( m - — 

lz zz V (J12 1)22 ■ 



(C(^7c - VJ - C(^jfc) + COjm)) ~ (COjfc ~ Vmk) ~ COjfc) + C(^mfc)) 
° m m (COtyfc ~ fijm) ~ CiVjk) + C{V>jm)) - (COjfc - A*mfc) ~ C(A»J*) + C(/W)) 

from which ([y]) follows. Substituting these results immediately yields fl42|) . 
After simplifying and again using fl3"3|) we obtain 

(83) 



Cjfc _ °"(^jfc) O'(fJ'jm) C(AW) ^(^jm + ^mfc — ^jfc) 



CjmCmk V\t*jk) V\yjm) ^Kymk) "Ky'jm ' H"mk H'jk) 

At this stage we have obtained 

l]x a(x - fi jm )a(x ~ Vjk + Hjr 



lM m (x) = -v 12 /(aOe^+fl"** 



a 2 (x)a(fi jm )a(fi jk - fi j7 



with a similar expression holding for l/ip m k and the desired expressions for A, r 
yi m fc and Ajfc. 

The final constraints arise by considering (^) which may be recast as 

G jm (x) - G mk [y) = din 1 - 



A jk (x + y) 

Now each side of this equation may be separately calculated and on comparison 
we find our last constraint. Using the form of Aj m , A mk and Aj k given by ( |42| ) 
we obtain for the left hand side 

C(X + fl mk - flj k ) - C{X - fij m ) - C(V + fijm ~ fljk) + C(V - Hmk), 

where substantial use has been made of ([34|). On the other hand, determining 
Gj m (x) directly from our expressions for ipj m yields 

G jm {x) = -Q{x - fj, jm ) - C(x + n jm - /M jk ) - F(x), 

where F(x) encodes the remaining functional dependence of ipj m . Comparison 
shows 

F(x) = -((x + fijm - Hjk) ~ ((x + ^mk - fijk) + const 

and so we have the final relation (E3f) 



G jm (x) =-£(x- fljm) +((X + flmk ~ fijk) + COnst. 
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We may for example choose 

<T(Hjk - ^jm)<?(Hjk - Aw)a 2 (x) 



/(*) 



a(x - jijk + Hj m )cr(x - n jk + \x mk ) ' 



Still f(x) is only determined up to a constant multiple of an exponential. This 
gives 

a{fi jm )a(x +fM- fx jm ) 

where we have set p, = —/J,j k + Hj m + Hmk an d k is an arbitrary constant. A 
similar expression holds for l/ijj mk . 

A. 2 Consistency for ( |25f) 

We next consider the constraints arising from the consistency of (22) when again 
Ajk, Aj m and A mk are given by (^) but now Aj m (x) = c 2 A mk (x). The latter 
of course means u jm = u mk , fi jm = fi mk , \ jm = A mfc and c jm = c 2 c mk . We now 
must find a function f(x), matrices U, V and A', A" such that 

c 2 A mk (x) \ f( x ) e -x*u ( $ ( x > v i* 



C m k(x) J \$'(x,U jk ) ' 

G m k{x)\ _ ,, x \"xy ( ®{x,Hjk) 

1 J \&{x,v jk ) 

We proceed in much the same manner as in the previous case and accordingly 
we will be less detailed. Now 

C2A mk {x) _ {x , +xll)x u n ®(x; Vjk) + u 12 $'(x; u jk ) 

1 v 21 ®(x; fi jk ) + v 22 ®'{x; fj, jk ) 

u 12 $(x; u jk ) [u n /ui2 + C(^jfc) - ({%) ~ ({vjk - x)} 
v 22 $(x; Ujk) [v 2 i/v 22 + C(n jk ) - C{x) - C(A»j* - x )] 

Again a comparison of zeros and poles shows 

U U /U 12 = ((Uj k - U mk ) - ((Ujk) + C(^mfc) 

V 2 i/v 22 = ((fi jk - Umk) ~ Citijk) + (iVrnk) 

Ujk ^mk f^jk l^mk 

A m fc + Cl^mfc) — Ci^mk) = Ajfc + ({Ujk) — Cil^jk) 

C 2 C mk V 22 = U 12 

V 2 \ U 21 



AjkX 



^mk ^jk 



V 22 U 22 
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These are the exact analogues of our earlier equations. In addition we have 
using the definition of C that 

C m k _ M2i$Q; Vjk) + u 22 <§'{x\ Vjk) __ A' mk _ 
A mk 2 u 11 $(x;v jk ) +u 1 2$'{x;v jk ) A rnk 

as well as 

Gmk __ VuMx] Mjfc) + V\2^'{X] fljk) 

1 v 2 i§(x; fi jk ) + v 22 §'(x; /j, jk ) ' 

Upon using substituting the form of A mk theses equations may be rearranged to 
give 

«2i + u 22 [((vjk) ~ C(x) - Ciyjk - x)} vn + vuiCifJ-jk) ~ C{x) - ({fijk - x)] 
^ 2 u u + u 12 [((is jk ) - C(x) - Q{v jk - x)\ v 21 + v 22 [C(fi jk ) - CO) - C(ji jk - x)} 

__ (T(x)cr(fJ, mk - l>mk)(?{.Hmk + "mk ~ X) 

a{u mk )a(x - v mk )o{ii mk )o{ii mk - x) ' 

As x — > we see 

u 22 v 12 _ 

C 2 1 A mk — U. 

«12 ^22 

Upon making use of this and further simplifying fl8~J|) yields 



^22,^21 , A/ _. N A/ _. N ,.,_. ^ <?[y mk )cr{y jk -v mk )(T{x-v jk ] 



C 2 ( + C(Wj-fc) - ((Vmk) - COjfc - *Vnfc)) , . , 

u 12 u 22 cr(Vj k )cr{x - v mk )a[x + v mk - u jk 

Vn.Vn o{\i mk )(j{iij k - n mk )a(x - fi jk ) 



^22 v 12 ■'■" ■ '' " <r(jtjk)<r(x — fJmik)<r(x + fin* — fijk) 

<j(v mk )cr(x - Vmk)<j(Hmk)cr(x ~ Hmk) 

As x — > /ijfc we find 

M22/M21 «ll^ 2 Kfc)ff 2 fe-M . 

c 2 — i ) — t — ^—77^ r ~~ ~ 

U 12 U 22 U 12 <j{Vj k )a{2v mk - U jk ) 

and asi-> z/jfc we obtain 

^12/^11 _ V21 , <y 2 {^mk)<y 2 {^jk - Vmk) _ , 

v 22 v 12 v 22 a(fi jk )a(2fi mk - n jk ) 

No further constraints are imposed on the matrices U and V. The ratio of these 
last two equations (taking into account that det U = Cj k det V) then yields 

Cjk _ v(Vjk) C 2 (/X mfc ) Cr(2u mk ~ Vjk) 

C 2 Cmk a (^jk) Cf 2 {Vmk) o{2\l mk - jX jk ) 

which is just the Vj m — > v mk limit of our previous result and we again obtain the 
relations stated in the theorem by similar analysis. 
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B Proof of Lemma 4 

Differentiating (|60|) with respect to y yields 

F'(x + y)= ( p(x)ij"(y)-(p'(x)ij'(y), 
and so upon letting y = we must solve the two equations 

F{x) = <j>(x)iJ/(0) - <p\x)ip(0), (85) 

F'(x) = (J)(x)ip"(0) - (f)'(x)ip'(0). (86) 

Now either F(x) ^ or F(x) = and 0(x) = ci exp(Ax), and if)(x) = c 2 exp(Ax). 
In the former case (perhaps by translating y if necessary) either i/)(0) ^ or 
■0'(O) 7^ and we consider these two cases separately 

If V'(O) y^ then utilising the group of symmetries of the functional equation 
we may set ^(O) = 1 and ^'(O) = 0. From (|55]) and (E3) 



F(x) = -<j)'{x) F'(x) = V"(O)0(x) 

and so 

<f>"( x ) = A 2 0(x) 

where A 2 = — ?/>"(0). Thus 0(x) = ciexp(Ax) + C2exp(— Ax) and F(x) = 
— A(ciexp(Ax) — C2exp(— Ax)). Therefore 

F(x + y) = — Xci exp X(x + y) + Xc 2 exp — X(x + y) 

= c 1 exp(Xx)(i/j'(y) - Xip(y)) + c 2 exp(-Xx)(ip'(y) + Xi/j(y)) 

which upon rearranging yields 

= ciexp(Ax)(^ / (y)-A^(y)+Aexp(Ay))+c 2 exp(-Ax)(^ / (y)+A^(y)-Aexp(A|/)). 

Thus either c\ = or c 2 = 0. The latter leads to 

= il>'(y)-\il>(y)+\exp(\y) 

and consequently (with the chosen initial conditions) 

^(y) = (l-Xy)exp(Xy). (87) 

which is of the required form. 

Finally, if ip(0) = 0, we may set ip'(0) = 1 and ^"(0) = 0. Now we have 
F(x) = 0(x) and F'(x) = — 0'(x), whence 0(x) = c\ is a constant. In this case 

tP(y)=c 2 + y, (88) 

again of the required form. 

□ 
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